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Hydrophobic hydration plays a crucial role in self-assembly processes over multiple length-scales, 
from the microscopic origins of inert gas solubility in water, to the mesoscopic organization of 
proteins and surfactant structures, to macroscopic phase separation. Many theoretical studies 
focus on the molecularly detailed interactions between oil and water, but the extrapolation of 
molecular-scale models to larger length-scale hydration phenomena is sometimes not warranted. 
Scaled-particle theories are based upon an interpolative view of that microscopic— >macroscopic 
issue. We revisit the scaled-particle theory proposed thirty years ago by Stillinger (J. Soln. Chem. 
2, 141-158 (1973)), adopt a practical generalization, and consider the implications for hydrophobic 
hydration in light of our current understanding. The generalization is based upon identifying a 
molecular length, implicit in previous applications of scaled-particle models, that provides an 
effective radius for joining microscopic and macroscopic descriptions. We demonstrate that the 
generalized theory correctly reproduces many of the anomalous thermodynamic properties of 
hydrophobic hydration for molecularly sized solutes, including solubility minima and entropy 
convergence, successfully interpolates between the microscopic and macroscopic extremes, and 
provides new insights into the underlying molecular mechanisms. The model considered here 
serves as a reference for theories that bridge microscopic and macroscopic hydrophobic effects. 
The results are discussed in terms of length-scales associated with component phenomena; in 
particular we first discuss the micro-macroscopic joining radius identified by the theory, then 
we discuss in turn the Tolman length that leads to an analogous length describing curvature 
corrections of a surface area model of hydrophobic hydration free energies, and the length-scales 
on which entropy convergence of hydration free energies are expected. 
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A. An empirical length characteristic of a liquid in 
coexistence with a dilute vapor 

The calling-out of a particular length-scale, and the 
role it plays in statistical theories of solutions, has been 
a prim ary feature of recent discuss ions of hydrophobic 
effects {Leel Il985t iLum et all Tl99^ . though the length- 
scales noted in those two cases were different from each 
other. It is interesting, therefore, to consider length-scale 
issues more b roadly. We can start by noti ng the classic 
suggestion of ijEgelstaff and Widoml Il970(l for a length- 
scale characteristic of a liquid in coexistence with a dilute 
vapor phase, namely the product of the liquid-vapor in- 
terfacial tension, 7, and the isothermal compressibility of 

the liquid, n T = — p. (^7) j see Fig. ^ The argument 
supporting this suggestion was physical, heuristic, and 
our discussion below of the scaled-particle theories will 
shed some additional light on this length-scale. It was 
immediately observed, however, that away from a criti- 
cal region jkt exhibited limited variation from liquid to 
liquid though the individual factors could differ by more 
than two orders of magnitude. 

Fig. H shows that the temperature dependence of the 
length jkt is qualitatively different for liquid water than 
it is for the organic solvents shown. This qualitatively 
different behavior is mostly ascribable to the fact that 
the compressibility of liquid water displays a minimum at 
46°C; it decreases with increasing temperature for tem- 
peratures lower than this, and has a smaller net variation 
over this temperature domain compared to the other sol- 
vents. Though liquid water is much less compressible 
than are the other liquids, the liquid-vapor surface ten- 
sion is much higher for water than for the other cases. 
This exemplifies the point that the product jkt has 
smaller variations than do the individual factors. Still, 
liquid water is distinguished from liquids generically by 
peculiar temperature dependences. Similarly, the pecu- 
liarities of hydrophobic effects are temperature depen- 
dences that imply entropic stabilization of conformations 
and assemblies of hydrophobic solutes in aqueous solu- 
tions. 

The theoretical description of hydrophobic effects has 
recently progressed markedly, and understanding the en- 
tropic interactions stabilizing of micelles, membranes, 
soluble proteins, and hierarchical biomolecular aggrega- 
tion in aqueous solution has similarly advanced. It is now 
recognized that scaled-particle theories can properly de- 
scribe primitive hydrophobic effects associated with the 
hydration of simple mono-atomic species. Scaled-particle 
theories identify — tentatively at first, but firmly as 
further information accumulates — a length separating 
microscopic from macroscopic descriptions of hydration 
structure. This establishes a radius at which microscopic 
and macroscopic descriptions of hydration structure can 
be effectively joined. Together with primitive constitu- 
tive information specific to liquid water, recognition of 
this joining radius provides an effective description of hy- 
drophobic effects for meso-scale aqueous structures. 



This paper traces those advances, specifically by lay- 
ing out the basic view generalizing applications of scaled- 
particle approaches. We establish the micro-macroscopic 
joining length, discuss the length — analogous to the 
Tolman length — associated with curvature corrections 
of a surface area model of hydrophobic hydration free 
energies, and finally examine the length-scales on which 
entropy convergence of hydration free energies are ex- 
pected. 



B. Sketch of experimental characteristics of hydrophobic 
hydration 

Hydrophobic phenomena usually do not occur in isola- 
tion from other interaction effects. The solutes that mo- 
tivate study of hydrophobic effects are typically molec- 
ularly complicated, water-soluble, amphiphilic chain 
molecules. 

Researchers studying these systems have been com- 
fortable, however, with a hydrophilic-hydrophobic di- 
chotomy. It is common to identify contributions 
to the hydration free energy above and beyond ob- 
yious hydrop hilic interactions as hydrophobic effects 
()Prattl Il998|) . This is particularly true if the tem- 
perature dependences of the comple mentary hyd r opho- 
bic interactions are also d istinctive l)Spolar et all Il989t 
ISpolar and RecordL Il994j) . A helpful review of hy- 
drophilic electrostatic interactions involved in protein 
molecular structure with an emphasis on the multi- 
ple length-scales involved appeared recently ijSimonsonl 
l2003j) . The present discussion emphasizes model solutes, 
inert excluded volume models, that permit study of hy- 
drophobic effects exclusively. Gases that are sparingly 
soluble in water, and small hydrocarbon molecules, can 
be brought into correspondence with hardcore molecular 
models. These models permit precision in isolating the 
temperature signatures that are the target of studies of 
hydrophobic effects. 

High precision calorimetric studies show that unfavor- 
able entropy changes upon transferring nonpolar solutes 
into water dominate hydration free energies at room tem- 
perature, and are only partly compensated by favorable 
enthalpy changes. These measurements are incongruous 
with regular solution notions, which ascribe insolubility 
to unfavorable cohesive interactions o f the solu t e wit h 
water compared to water with itself ijLazaridisl |2001) . 
Of foremost importance, the entropies and enthalpies 
of hydrophobic hydration are strongly temperature de- 
pendent, following from the large positive heat capac- 
ity changes observed. Thus, at elevated temperatures 
the roles of entropy and enthalpy are reversed, with un- 
favorable enthalpies dominating hydration free energies, 
partly compensated by favorable entropies. The result- 
ing solubilities of nonpolar gases are nonmonotonic, ex- 
hibiting a solubility minimum between T = 310 K and 
350 K. Ana logously, prote i ns undergo hot and cold de- 
naturation ((Brandts!. Il964t iFranks and Hatlevl Il99l|) as 
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noted above, while ionic and nonionic surfactants display 
a minimum in their crit ical micelle c oncentra tions with 
respect to temperature l)Chen et all Il998albl ) , pointing 
to a common underlying mechanism with the solubility 
behavior of nonpolar species. 

The experimental entropy changes for a range 
of hydrophobic solutes inters ect at v alues close 
to z e ro near T = 385 K llBaldwinL Il986t iLee 
19911 iMakhatadze and Privalovt Il995t iMulleit 1199.1 



Murphv et all 1199(1 iPrivalovt Il979t iPrivalov and fid 



1988). The coincidence of this entropy convergence tem- 
perature for hydrocarbons with comparable behavior in 
protein unfolding has been used as empirical justifica- 
tion for the hydrophobic core model for protein folding, 
and has influenced the interpretation of biomolecular as- 
sembly. But in the complex context of soluble protein 
molecules, a cl ear relevance of entropy conve rgence can 
be questioned i|Robertson and Murnhvl Il997h . 

It has been traditionally argued that these hydropho- 
bic entropy effects stem from orientational constraints 
on water molecules in the hydration shell of nonpo- 
lar solutes, constraints that maintain the integrity of a 
hydrogen-bonding netw ork forming cage-l i ke str uctures 
or microscopic icebergs ijFrank and EvansL Il945[l . The- 
oretical studies of the impact of local clathrate for- 
mation about krypton in water conclude that literal 
clathrate structures are far too unf avorable to play 
a role in hydration thermodynamics l|Ashbaugh et all 
2003). Experimental probes of the local structure of 
water proximal to purely nonpolar solutes are scarce, 
and hampered by the low solute-concentrations attain- 
able. The structures that have been studied by neu- 
tron and X-ray scattering techniques suggest that while 
water surely adopts orientational preferences in the hy- 
dration shell of nonpolar moieties, the solute-induced 
structu re is more diso r dered than in ice or clathrate hy- 
dratesl|BcjwrQri_g^a/] 1 ^ Neilsonl 
ll994HDe.Tong et n,l\ . ll 997HFilir)nom et oil ll 997ft . 

Molecular-level investigations of hydrophobic hydra- 
tion have been la rgely theoretical and simulation cf- 

onrv 



forts llHendersonl l2003 iHummer et all Il998l [l99fi : 
Pierottil Il976t iPohorille and Prattl Il990t iPrattt 12002 ; 



Pratt and PohorilleLll992tlStillingeIll973j) . Water struc 
ture in the vicinity of hydrophobic species, includ- 
ing orientational preferences, has been connected to 
the entropies of hydrop hobic hydration by a correla- 
tion function expansion ll Ashbaugh and Paulaitisl Il996t 
lLazaridis and Paulaitisl Il992t ISilverstein et all l200l[) 
Pratt and Pohorille, on the other hand, demon- 
strated that the low solubility of atomic solutes in 
water arises from the narr ow distribution of cavity- 
open ing fluctuations in water l|Pratt and PohorilleLll992l 
Il993l). With connections t o the Pratt-Chandler theory 
l|Pratt and Chandlerl Jl977j) and its Gaussian field inter- 
pretation ijChandleit Il993t iLvrm^^^ 1999) , an infor- 
mation theoretic (IT) model ijHummer et a/.lll998ll2000|) 
then provided a quantitative link between the micro- 
scopic density fluctuations determined from water oxygen 



pair correlations and the hydration free energies of hard 
solutes. More importantly, the IT model implicated the 
unusual equation-of-state properties of water as a domi- 
nant factor in hydrophobic hydration, differentiating wa- 
ter from other common solvents. 

In addition to capturing temperature and pressure ef- 
fects as sociated with hydroph obic hydration and inter- 
actions ijHummer et all Il998|) . the IT model provides a 
facile description of the enhanced solubility of nonpolar 
species in D2O compared to H2O as a result of differ- 
ences in the isothermal co mpressibilities of these isotopic 
alternative forms of water <)Hummer et a/.U2000() . 

The latter point suggests again the length ^kt of 
Fig. ^ because of the involvement of the measured 
isothermal compressibility with its specific temperature 
dependence. But the product involves the surface 
tension 7 also, and the surface tension has been invoked 
in empirical correlations connecting measurable proper- 
ties of liq uid water and h ydrophobic effects. Years ago, 
however, l|Tanfordl Il979j) pointed out the large discrep- 
ancy between the measured water-hydrocarbon interfa- 
cial tension and the effective microscopic surface tensions 
obtained from hydrocarbon solubility data. A correspon- 
dence between macroscopic and microscopic surface ten- 
sions has been contentious because of their fundamen- 
tally different temperature dependence. As noted above, 
it is temperature dependences that have been the basic 
puzzles of hydrophobic effects. 

More recent examples of the distinction between 
molecular and macroscopic hydrophobic interac- 
tions are found in measurements of the long- 
range attractive fo rce between macroscop ic hy- 
drophobic surfaces llChristenson and Claessonl Jl988t 
llsraelachvili and Pashlevl Il982t iPashlev et all Il985|) 
which have not been explained on the basis of molec- 
ular hydrophobic effects. Vibrational sum frequency 
spectroscopy suggests that hydrogen-bonding of water 
molecules is weaker at macroscopic water-carbon tetra- 
chloride and water-hexane interfaces than near individual 
hydro phobic species dissolved in water ijScatena et all 
I2001D . The lack of a definitive interpretation of these 
surface force measurements, and the changes in water 
energetics at macroscopic interfaces underscores the need 
for a quantitative theory beyond molecular hydrophobic 
effects. In general, the need for a unified, quantita- 
tive description of both molecular and macroscopic 
hydrophobic phenomena arises because hydrophobic 
driving forces play an important role in self-assembly 
on multiple length-scales and the fact that quantitative 
descriptions of these driving forces are derived from 
molecular solubility data, macroscopic interfacial tension 
o r interpolations of these quantities 
. Il999t [ Ashbaugh and Paulaitisl 120011: 
l2000t Irlermannl. Il977t ISharp et all 



measurements, 



(Ashbaugh et al 
G alli cchio et all 
Il99li lTanfordlll97: 

l)Lum et all Il999j) suggested bridging these disparate 
length-scales by incorporating a Gaussian field theory for 
molecular level fluctuations with mean-field theory for 
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larger scale structures ultimately responsible for macro- 
scopic phase transitions. Their approach successfully 
predicts many of the thermodynamic anomalies char- 
acteristic of small molecule hydration, and goes fur- 
ther, predicting the onset of long-range hydrophobic 
forces between surfaces as a result of an aqueous liquid- 
vapor phase transition in confined geometries. Indeed 
surface force apparat us studies of the long-range hy- 
drophobic interaction ijChristenson and Claessonl Il988j) 
and simulations of water confined between repulsive 
oblate elli psoids observed ca vitation between nonpolar 
surfaces ijHuang et all l2003|) . consistent with theoreti- 
cal predictions. Mean-field modeling and simulations of 
methane clusters, however, suggest that when ubiqui- 
tous attractive interactions between water and hydropho- 
bic surfaces are taken into account, surface and confine- 
ment induced local stru c tural changes are suppressed 
llAshbaugh and Paulaitisl l200lt lAshbaugh et all 12004 
[ChauL l2003tlDzubiella and HansenLl2003tlTruskett et all 
l200l|) . Moreover, experiments on the effects of electrolyte 
addition and degassing on the range of surface forces, and 
the stability of surfactant free a queous emulsions chal- 
lenge the theoretical predictions (IConsidine et all fl99£ ; 
Kokkoli and Zukoskilll99llPashlevll2003l:lPashlev et al . 

A conceptual basis for unifying molecular and macro- 
scopic hydrophobic hydration can be found in scaled- 
particle theory (SPT). Thirty years ago iStillingeilll97l 
presented an influential p aper on the applicat i on of 
the c l assic SPT o f Reiss ijPierottl Il976t iReissl Il965l 
Il977t iReiss et aR Il959j) to the hydration thermody- 
namics of purely excluded volume solutes. The pur- 
pose of that paper was, in part, to illuminate the 
pitfalls and difficulties in applying classic SPT, origi- 
nally developed for hard-sphere fluids, to aqueous sol- 
vents l)Ben-Naim and Friedmanl Il967f) . In doing so 
Stillinger opened up new avenues of inquiry into hy- 
drophobic hydration within the context of SPT. Nev- 
ertheless, direct exploration of the validity and con- 
sequences of Stillinger's revised theory have been rare 
l|Pratt and PohorilleL Il992l Il993ft . We presently revisit 
SPT and critically discuss its implications in light of our 
current understanding of hydrophobic hydration. For the 
first time, we demonstrate that the revised SPT repro- 
duces many of the characteristic thermodynamic signa- 
tures of molecular hydrophobic effects and can be used to 
extend the results of molecular simulations of small hard 
hydrophobic solutes in water to meso- and macroscopic 
surface hydration. The present analysis provides insights 
into the differences and similarities for hydrating molec- 
ular and macroscopic surfaces. In addition, we examine 
the validity of surface area correlations commonly used in 
biophysical models for hydration thermodynamics over a 
range of length-scales, as well as the origins of entropy 
convergence behavior at molecular lengths scales and how 
solute size moderates the convergence temperature. 



II. A PRIMER ON SCALED-PARTICLE THEORY 

A. Classic scaled-particle theory 

The chemical potential of a hydrated mono-atomic so- 
lute, 'A', can be expressed formally as 

PA (aq) = kT In [p A (aq) A A 3 ] + Ma (*0 (1) 

where we have adopted Ben-Nairn's standard 
state in the definition of the chemical potential 
l)Ben-Naim and Marcusl Il984j) . Here k is Boltzmann's 
constant, T is the temperature, pa (aq) is the solute 
number density in the solution, Aa is the thermal de 
Broglie wavelength of the solute, and /j, A x is the excess 
chemical potential, i.e., the coupling work of turning 
on interactions between the solute and water, which 
disappears if those interactions vanish. At thermal equi- 
librium the Ostwald partition coefficient determining 
the distribution of the solute between an aqueous and 
ideal gas \p, e * (ideal) = 0] phase at infinite dilution is 
(|PollackL fl9^ 

K cq = PA . (aq) = exp (aq) /kT] . (2) 

Pa (ideal) 

Thus, the excess chemical potential is central to resolv- 
ing the aqueous solubility of the solute. Confining our 
discussion to impenetrable hard sphere (HS) solutes, the 
solute excess chemical potential in water is 

MA x (aq) = -kTlnp (R) (3) 

where the insertion probability, po(R), is the probability 
that a solute-sized stencil randomly placed in water is 
devoid of water oxygen centers; this follows directly from 
Widom's potential distributio n theorem for species inter- 
actin g with a hard potential ijPratt et all Il999l IWidomL 
Il982t) . The solvent accessible radius, R, is the radius 
of closest approach between the solute center and a wa- 
ter oxygen. The solvent accessible radius is the sum of 
the van der Waals radius of a water molecule and the ra- 
dius of the hard-sphere solute, i.e., R = (er\vw + ct aa) /2. 
The insertion probability is the fractional free volume 
offered by the solution — or the available volume in 
a mo r e spe cialized language attributed to Boltzmann 

jm[l98i), i-e., Po(R) = VfceeAW 

Figure [3 gives a physical picture of this available vol- 
ume. Given a molecular snapshot of liquid water — the 
left box in Fig. [21 — Vf rec is determined by the points 
at which a hard sphere could be successfully implanted. 
Graphic display of those successful placements thus yields 
a negative image — the right box in Fig. [2] — of the 
molecular configuration from which the analysis started. 
Vfrec decreases with increasing cavity radius, as success- 
ful insertions become increasingly rare. po(R) is subse- 
quently determined as an ensemble average over a large 
number of molecular configurations. 

An alternative relationship for the chemical potential 
that draws a connection to the structure of water in the 
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vicinity of the solute is 



ff* (aq) = kT / o w G(r)4 7 rr 2 dr 
Jo 



(4) 



where pwG(r) is the density of water in contact with a 
hard solute of contact radius r (|StillingerL Il97^ . The 
product kTpwG(r) has units of force/area, and is the 
pressure due to solvent collisions with the hard surface of 
the solute. The chemical potential then is the reversible 
pressure- volume work required to expand the solute into 
the solution. The etymolo gy of scaled-partic le theory 
derives from this expression (Rciss, ll965Ul977h since the 
solute is introduced by scaling-up from a particle with 
a solvent accessible radius of zero to a final size of R. 
The relationship between the contact function and the 
insertion probability is 



G{R) 



dln Po (R) 



AttR 2 Pw OR 



(5) 



as determined by differentiation of Eq. (0J. For suf- 
ficiently small solute cavities that only one solvent 
molecule might fit within the cavity boundary, the in- 
sertion probability is 



47T 

p (R^Q) = l-—p w R 3 . 



(6) 



This is merely Vf ree /Vtotai m the case that the solute is so 
small that it never interacts with more than one solvent 
molecule exclusion sphere at once. Then the excluded 
volume is (47r_R 3 /3) Nyj-; the more general case comes- 
up below. The expression corresponding to Eq.|Hlfor the 
contact function is 



1 



(7) 



These expressions are accurate in solvents with realistic 
interactions for cavities up to a radius of R rs ctww/2, 
and are exact in hard sphere fluids then because the so- 
lute corresponds to a point particle with a van der Waals 
radius of zero. Larger hard-sphere solutes can contact 



more than one solvent molecule at once, and addressing 
larger solutes requires information on dauntingly compli- 
cated multi-body correlations. 

In the limit of a macroscopically large cavity, the con- 
tact correlation function can be represented as an asymp- 
totic expansion in 1/R 



G(R) 



j>0 



R> 



(8) 



Retaining contributions up to j = 2 yields an expres- 
sion equivalent in form to that required by classical 
ther modynamics for the force a cting on the cavit y sur- 
face (iHendersonl 120021: Ipierottil 119761 iReissl I1965L Il977t 
IStillingerL ll97^T 



kT Pw G{R) ~ p 



27oo 47^(5 



R 



R 2 



(9) 



where p is here to equal the liquid saturation pressure 
Psati 7oo is henceforth the surface tension of a flat in- 
terfac e that was denoted by 7 above, and <5 ijTolmanl 
1949) describes the initial curvature correction to the 
surface tension. a The coefficient G3 is zero so that 
the c hemical potential is free of logarithmic contribu- 
tions (Stil linger and CotterLll97lh though analyses of the 
possibility of logari t hmic corrections are still of interest 
IjEvans et all 120041 12003(1 . Higher order terms in the 
asymptotic expansion are not generally ava ilable. Such 
considerations motivated l|Reiss et all [19591 to truncate 
Eq. (JHJ after the initial curvature correction in order to 
develop a tractable, physically reasonable model for the 
contact function. 

a The notation of Eq. follows a typographic confusion wide- 
spread across the present problem. S here is one-half the conven- 
tional Tolman length. Expressed more basically th e final term 
on Eq. 151 is f27oo /R) times the Tolman length. See iHendersonl 
l2002t iMoodv and Attardl |2001|) . We thank J. R. Henderson for 
a discussion of this point. 

Evaluating the j = term in Eq. © with the mea- 
sured equation of state, and the j = 1, 2 terms in the 
expansion by requiring the microscopic and macroscopic 
limits smoothly meet at <7ww/2, yields 



G(R) 



Psat 

kTpw 



i-4rPw-R 3 

2+1) 2p B . 



kTpv 



( £ww\ 
\ 2R J 
2 



(1+21?) 



Psat 



fcTpw 



( <7WW 

V 2R 



R < (Tyyw/2 



R > <Xvyyv/2 



(10) 



where r\ = ^pw^ww 3 is the solvent packing fraction. Integration of the contact function yields the excess chemical 
potential 



Infl 



r -PwR 3 



Ma 
kT 



■ln(l-r?) + 



2(1-7,)" 



kTpv 



3r,(2+r,) _ 3r,p 3a 
2(l-r,) 2 kT Pv 



2R 



3r,(l+2r,) 



kT pw V cww 



kTpv 
2R 



R < cr ww / 2 
R > <7ww/2 ■ 



(11) 
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Eqs. (|10(l and constitute the cl assic SPT or i ginall y 
developed for the hard sphere fluid (jReissL Il965l Il977l). 
but which was s ubsequ ently applied to water bv (|Pierotti| . 
1 «»—(»:. and |Lei[l985). Fundamental difficulties arise in 
the application of classic SPT to water, however, includ- 
ing the erroneous prediction that the surface tension of 
water increases with tempe rature and pass es through a 
maximum near T — 425 K (Stillinacr, 1973). 



Waals diameter, ctwWj that might differentiate water 
from other solvents, and thereby limits the interpreta- 
tion of complex hydration phenomena. To consider this 
problem more generally, we note that the insertion prob- 
ability is formally expressed on the basis of s olvent struc- 
ture by an inclusion-e xclusion development ijReissi 1 19651 
119771 IStillingerl 11973ft 



B. Revised scaled-particle theory 

The scaled-particle model described above incorpo- 
rates little molecular detail beyond the assigned van der 



Po (R) = l + J2 { -^ [ ■■■[ S W(n...r n )d 3 r 1 ...d 3 r tt (12) 

where V(R) = AttR 3 /3 is the observation volume, and (ri . . . r„) are the n-body solvent oxygen distribution 
functions. The terms in this series vanish for n exceeding the maximum number of solvent molecular centers that 
can be packed into a sphere of volume AttR 3 /3. It is on this basis that the limiting result, Eq. is established. 
As noted above, these distribution functions are complicated, not routinely available beyond t he pair-distribution 
function (n = 2), and there has been only one detailed investigation of terms beyond 2nd order ijGomez et all 119991 
iPratt et aZJ . Il9 99) . Considering the small cavity pair-correlation (n — 2) contribution and the as ymptotic macros copic 
thermodynamic limits, Stillinger proposed a revised expression for the cavity contact function l)Stilringerlll973[) 

l+^Lti R 9 {2) (r)r 2 (r-2R)dr R < R* 

G(R) = { l-4fpwi? 3 + (^)"/ 2i? S< 2 )W(r3/6-2^,. +8fl 3 /3))dr > - ( 13 ) 

Peat I 27QQ i'f^S I X R ^> R* 

fcTp„ T kTp w R kTp w R 2 T R 4 - ' Jt ^ " 



where R* is the radius at which n = 3 correlations be- 
gin to contribute to the cavity insertion probability of 
Eq. 1121 While the experimental pressure, surface ten- 
sion, density, and solvent radial distribution function are 
employed, 5 and A are treated as adjustable parameters 
chosen so that the small cavity and macroscopic limits 
of the contact function join smoothly at R* . This ex- 
pression incorporates molecular information on the pair 
structure of water as well as the known macroscopic 
properties of bulk water and its interfacial behavior, and 
thereby is ex pected to discriminate more sen sitively be- 
tween w ater jAshbau gh and Patilaitis , l200l[) and other 
solvents ijHuang and Chandlerl 120006 ). Indeed, it has 
been demonstrated by extensive molecular simulations 
that Eq. (|13fl provides a description superior to the clas- 
sic SPT expression Eq. (|10fl of the solvent contact density 
for solutes several times larger than the solvent. 

Stillinger's revised SPT prediction for G(R) relies 
on the assumption that multi-body water correlations 
at intermediate, but molecule-sized, solute radii are 
adequately represented by the parameters 5 and A 
fitted at a radius R* . Numerical experimentation 
ijPratt and PohorilleL Il992[) with this parameterization 



shows that the revised SPT is most sensitive to the pa- 
rameter R* . This parameterization might be improved 
by involving results over a range of radii, including so- 
lute sizes for which multi-body correlations are signif- 
icant. While this information is not readily available 
experimentally, multi-body correlation contributions to 
the hard-sphere solute chemical potential can be interro- 
gated by direct evaluation of the insertion probabilities 
from molecular simulations of water. In the spirit of Still- 
inger's revised SPT, we interpolate between the chemical 
potential evaluated for molecular length-scales from sim- 
ulation, and the macroscopic thermodynamic limit 

A4 x (i?) = -kTln Po (R)\ sim f(R) 

+ Ma 0R)Ucro[l -/(#)] (14) 

where f(R) is a switching function equal to one below 
i? s i m and zero above i? ma cro, smoothly interpolating be- 
tween these between these two limits. Presently we use 
a cubic spline interpolating function, though other rea- 
sonable functions yield essentially indistinguishable pre- 
dictions. The macroscopic chemical potential, deter- 
mined by integration of the macroscopic cavity expres- 
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Eq. 



III. APPLICATION TO HYDROPHOBIC HYDRATION 



R 



e — 167ri?7 oc (5 



4irR 2 loc + — R 3 Ps , t (15) 



Rather than fitting the microscopic and macroscopic lim- 
its at a single point as in Eq. (|13[) . the parameters S, A, 
and the integration constant e are fitted to the simulation 
results between i? s im and i? m acro- The contact correla- 
tion function is then determined by differentiation of the 
chemical potential [Eq. J5J]. An additional benefit of fit- 
ting Eq. (|15|l to the simulation insertion probabilities is 
that we do not have to evaluate numerically first and sec- 
ond derivatives of the simulated insertion probabilities, 
which become more statistically uncertain with increas- 
ing cavity size. 



C. Computational implementation 

The computational implementation follows standard 
procedures for sampling molecular configurations of 
liquid water, and builds f rom the original work of 
llPohorille and Prattt Il990t iPratt and PohorilleL Il992t 
1993) in evaluating cavity statistics therefrom. Water 
configurations were generated using Monte Carlo sim- 
ulations in the canonical ensemble ijFrenkel and Smil 
2002). Bulk water was modeled using 268 SPC/E 
water molecules with periodic boundary conditions 
iBerendsen et oil 119871) . SPC/E was chosen because it 
provides an accurate representations of the structure, 
equation-of-state, and interfacial tensi on of liquid wa- 
ter over a broad range o f temperatures ijAleiandre et all 
ll995UHura et nil 2003). Lennard- Jones potential inter- 
actions were evaluated smoothly truncating the potential 
based on the separation of water oxygen atoms between 
9.5 A and 10 A, while longer ranged electrostatic interac- 
tions were calculated using E wald summation with con - 
ducting boundary conditions l)Frenkel and Smitl |2002|) . 
Simulations were carried out from T — 260 K to 470 K in 
10 K increments at the experimental liquid density along 
the saturation curve and into the supercooled regime 
l|Hare and Sorensenl Il986ll . After an equilibration phase 
of at least 10 5 MC passes (where one pass corresponds to 
one attempted move per water molecule with 30% move 
acceptance), 5 x 10 6 MC production passes were carried 
out for analysis of thermodynamic averages. After each 
50 MC passes, 10 5 particle insertions were attempted to 
estimate po(R), so a total of 10 10 insertions were at- 
tempted at each temperature. Statistical uncertainties 
were determined by grouping results into block averages 
over 10 6 MC passes each. 



A. Cavity contact functions: the micro-macro joining 
boundary 

The cavity contact function at T = 300 K is shown 
in Fig. |3J Beginning at a value of one at zero radius, 
the cavity contact density increases with increasing R. 
Simulation values of G(R) appear to plateau at a max- 
imum near 3 A. Just beyond this radius, the simula- 
tion results for G(R) become progressively noisy as a 
result of poor sampling of infrequent large cavity fluc- 
tuations. Detailed calculations for specific values of R 
greater than 3 A ijAshbaugh and Paulaitisl l200l|) have 
established that this is indeed the region of a maximum 
in G(R), and that G(R) is qualitatively described by 
Stillinger's revised scale particle model. A dominating 
observation is that this curve imposes a non-arbitrary 
definition of a length-scale for the present problem: the 
radius i? ma x at which G(R) is a maximum. Solutes with 
smaller radii are identified as intrinsically microscopic in 
scale. The description of larger solutes can be built from 
a macroscopic perspective. An interpolative strategy ex- 
tending to large solutes, such as that adopted here, is 
likely to be effective if the region at which the molecu- 
lar and macroscopic expressions are joined encompasses 
-Rmax- The revised SPT fit, determined by differentiation 
of Eq. l(T4*|) fitted to the simulation insertion results be- 
tween i? s j m = 2.5 A and i? maC ro = 3.5 A, extends G(R) to 
R larger than observed directly. The revised SPT result 
places the maximum contact density at R = 3.0 A where 
G(R) w 2.3. Solutes of this size are candidates for most 
hydrophobic because the compressive pressure exerted by 
the solvent is largest in this case; see Fig. 0] 

Beyond this maximum, water pulls away from the cav- 
ity surface with increasing size. At the size R « 10 A the 
contact density equals the bulk density of water, decreas- 
ing further for larger cavities. In the limit R — > oo, the 
contact correlation approaches 



Psat 

p w fcT 



2 x 10 5 for wa- 



ter at T — 300 K, and the pressures here are sufficiently 
low that they do not influence the contact function for 
molecular and mesoscopic cavities at any of the temper- 
atures considered. 

This appar ently anomalous drying behavior was an- 
ticipated by (jStillingeii Il973j) . and has only recently 
been confirmed by molecular simulations in Lennard- 
Jones and aqueous solvents fAshbaugh and Paulaitisl 
l200li iHuang and Chandleil l2000a|) . The surface drying 
has been previously interpreted in terms of an effective 
exp ulsion potential between water and the solute cav- 
ity Jlum^eLanTGardi GUI E^ekTeliil EI!!)- In 
bulk water, the individual water molecules feel attrac- 
tive interactions with the other water molecules, and the 
average force on a water molecule in the bulk is zero. To 
approach a large solute, however, a water molecule must 
shed hydration partners and/or limit their possible con- 
figurations. This unbalances the interactions with the 
aqueous medium and gives rise to an additional repul- 
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sive force between a water molecule and the surface. If 
the solute is unable to compensate for these lost interac- 
tions to counter the cavity expulsion potential, water is 
repelled by the surface. 

The classic SPT prediction for the contact function 
[Eq. HI L)|l ] is in qualitative agreement with the simulation 
and revised SPT results [Fig.|3j. Classic SPT predicts a 
maximum in the contact density, followed by a decrease 
to values below the bulk density of water with increas- 
ing cavity radius. The quantitative agreement is poor, 
however, even if the effective diameter, owwi of water 
is treated as an adjustable parameter. Notably, clas- 
sic SPT predicts that the maximum in G(R) is shifted 
to smaller radii of R ~ 2 A. Water has an open struc- 
ture favoring larger cavities at these packing fraction s 
ijPohorille and Prattt Il990t iPratt and Pohorillel Il992|) . 
The resulting maximum in the pressure acting on the 
solute surface, kTpyjG(R), then is shifted out to larger 
cavity radii for water. 

If the objective of classic SPT is to reproduce the chem- 
ical potentials of solutes using Eq. (f 1 1|1 up to radii of R « 
3.3 A, encompassing the sizes of a number of nonpolar 
gases, a typical diameter , trww, of water fitted to ex- 
perimental data is 2.7 A l)LeeUl985j) . though a more ap- 
propriate value based on the simulation results reported 
in Fig. |3 is 2.8 A. In this case, the fitted radius of wa- 
ter splits the difference between the over-prediction of 
G(R) at intermediate radii by the classic SPT, and the 
under-prediction of G(R) at radii close to the maximum 
solute size to balance out these inaccuracies in the calcu- 
lation of the chemical potential. This oww assigned to 
water molecules then is a consequence of the fitting, and 
does not contribute to the interpretation of the molec- 
ular signatures of hydrophobicity. Indeed, if we extend 
the predictions of classic SPT outside the range fitted 
for small solutes, we find the theory under-predicts the 
hydration free energies of mesoscopic cavities, and that 
a larger water diameter — ctww ~ 2.9 A! — is required 
to match the drying observed in G(R) [Fig. . 

While the variation of this size parameter may seem 
small, the chemical potential depends on the integral of 
Eq.0] and small differences in crww significantly alter the 
predictions. Broadly viewed, this is the natural observa- 
tion that slight adjustment of boundary information in 
boundary value problems can make large changes away 
from the boundary. 

Revised SPT predictions for G(R) as a function of tem- 
perature are shown in Fig. |SJ While all the curves are 
qualitatively similar, the magnitude of the contact func- 
tion decreases with increasing temperature. Classic SPT 
fails to describe this temperature dependence of G(R). 
Fig. [5] also indicates that the length defined by the maxi- 
mum of this curve decreases with increasing temperature, 
consistent with the suggestion of Fig.^ In following sec- 
tions we use these revised SPT results to draw conclu- 
sions regarding the size and temperature dependence of 
the hydration of hydrophobic hard spheres. 

The discussion in Sec. II. Al of the length-scale joo kt 



sheds some light on the results of Fig. [3] Specifically, we 
can put together the crudest of models of the hydration 
free energies for the small-scale and large-scale problems, 
and in that way get a crude characterization of the length 
at which these different descriptions for G (A) match. For 
the small-scale — but not the smallest-scal e as in Eq.Efl 
the information th eory approach suggested l|Garde et all 
ll99fitlPrattl 120021) (3^ « (n) Q 2 /2 (8n 2 ) Q . If we evaluate 

everything on a macroscopic basis, then G (R) w 2pwK T 
from Eq.[S]. This is an extremely crude estimate, but we 
rely only the trends with changes in temperature in the 
present argument. For the large-scale problem, we use 
Pn e j£ ~ 47ri? 2 /?7 00 , expecting that the pressure will be 
negligibly low for this consideration. Then G (R) ~ jj 2 -- 
These two estimates match when i? max ~ ^loo^T- The 
first of these estimates, that from the information theory 
contribution evaluated on macroscopic information, is 
too large by about a factor of 4. Exploiting an empirical 
factor of 4 would bring the estimate i? m ax ~ 4 x 2700 kt 
into a realistic range of molecular sizes. But the present 
argument is crude enough that only the dimensionally 
natural linear correlation with the length 700 «t should 
be taken seriously. 



B. Physical relevance of hard-core model solutes to 
structural theories of hydrophobic effects 

Hard core models of solute-water interactions stud- 
ied in this context, which are the b a sis of scaled- 
particle theories, are motivated ijPrattl l2002f) by the 
broa d success of van der Waals theories of liq- 
uids llBarker and Hendersoiill97fiUChandler et oJ.l . ll983l 
iLebowitz and W aisman. ll980HWidomLll967|) . Physically 
expressed, van der Waals approaches are appropriate for 
the description of dense liquids because the disorder and 
high density can limit structural fluctuations to length- 
scales associated with the variation of intermolecular re- 
pulsive interactions, small compared to the range of inter- 
molecular attractive interactions. When attractive inter- 
molecular interactions are weak on a thermal scale, that 
helps too since a van der Waal approach treats those in- 
teractions perturbatively. Then an approach which sep- 
arates the effects of intermolecular repulsive and attrac- 
tive interactions can be successful. In that case, repul- 
sive interactions present the first challenge to theories, 
and hard-core interactions are natural simple models for 
those excluded volume interactions. This is the argument 
for the physical relevance of hard-core models solutes to 
theories of hydrophobic effects. 

With this background, the most important physical 
observation on the large R behavior of the G(R) results 
of Fig. is that those results would be expected to be 
sensitive to attractive solute-water interactions, if they 
were to be included. When R is large, the local density 
in the vicinity of the hard sphere solute can be low, and 
the argument above that fluctuations do not access the 
length-scales comparable to the range of natural attrac- 
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tive interactions does not apply. Simulation evidence 
does support the view that results can be sensitive to 
inclusion of natural attractive interactions when the so- 
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Nevertheless, hard-core models of solute- water interac- 
tions serve as an valuable reference point for at least two 
reasons. A first reason is conceptual and reductionist. 
This simplified case has historically been considered as 
expressing the basic puzzle of hydrophobic effects. (The 
extent to which that is true is one of the issues addressed 
here!) Solving this basic puzzle enables specific cases to 
be described by combination of what is understood for 
the simpler cases. A second reason that hard-core models 
of solute- water interactions are valuable is that for R not 
too large the results should be less sensitive specifically 
to the case of physical interest. (Support for this view is 
noted at the appropriate places in the succeeding discus- 
sions.) From this point of view, then, the careful study 
of the large R behavior of hard sphere G(R) assists in 
refining the description of intermediate R behavior, in- 
cluding the region of the maximum corresponding to the 
most hydrophobic solutes. 



C. Hydration thermodynamics of hydrophobic species: 
temperature signatures and solubility minima 

The hydration free energy of a methane sized (R = 
3.3 A) hard-sphere solute in water as a function of tem- 
perature along the saturation curve is shown in Fig. [5] 
The simulation results for the chemical potential pass 
through a maximum at T w 400 K, at which point the 
hydration entropy defined by s A x = — <9^ A x /<9T| sat is zero. 
To extract the enthalpy and entropy of hydrophobic hy- 
dration from the chemical potential, we assume that the 



heat capacity dh c £/dT\ s 



Tds c l/dT\ s 



(T) is 



independent of temperature. In this case, the hydration 
enthalpy, entropy, and free energy are 

hf = h% (T ) + (T — T ) cf (T ) , (16a) 
*T = »? (To) + In (To) , (16b) 

= nf (To) + (T - T ) (4* (T ) - 4 X (To)) 

-rin^W(To), (16c) 

respectively, and ^ (T ) = h% (T ) - T sf (To). 

The enthalpy and entropy of hydration of the methane- 
sized hard-sphere solute solute is shown alongside the 
fitted free energy in Fig. HJ The hydration cn- 
tropy is negative and unfavorable at room temperature 
llBlokziil and Engbertslll993t lK auzman nl ll959tlTanfordL 
Il980|) . With increasing temperature the entropy changes 



sign indicative of a positive heat capacity increment. 
The entropy and heat capacity at T = 298 K for the 
hard-sphere solute are —69.5 J/(mol K) and 214 J/(mol 
K), respectively, which is in excellent agreement with 
the experimental values for the entropy and heat ca- 
pacity of -66.7 J/(mol K) and 209 to 237 J/(mol K) 
for methane at T = 298 K llBen-Naim and Marcusj . 
[ 19841: lLazaridis ; and Paulaitisl Il992l iNaehibi ef all 119861 
iRettich et all Il98lj) . Over most of the temperature 
range considered, the hydration enthalpy is positive and 
unfavorable for hydration, in disagreement with the ex- 
perimental enthalpy for methane of —11.5 kJ/mol at T = 
298 K, largely a result of the neglect of attractive inter- 
actions with^ater^^he iceftero hypothesis of Frank and 
Evans ijFrank and EvansLll945|) suggests that local freez- 
ing of water in the vicinity of hydrophobic species con- 
tributes to the experimental negative hydration enthalpy. 
In the case of the methane sized hard sphere though, the 
hydration enthalpy at T = 298 K is 5.0 kJ/mol, contrary 
to customary expectations. 

The ratio of the chemical potential and kT dictates 
the Ostwald solubility. — In K eq for the methane sized 
solute in water as a function of temperature is shown 
in Fig. \7\ This quantity passes through a maximum 
near T = 280 K, corresponding to a minimum in the 
solubility. Thi s observation is in agreement with infor- 
mation theory l| Garde_i st all. Il999fl and equation-of-state 
ijAshbaugh et a7] . l2002y ^nodels of hard-sphere solubili- 
ties that link the solubility minimum to the density max- 
imum at T — 277 K for pure water. The solubility min- 
imum corresponds to the point at which the enthalpy, 
^cx _ —T 2 d(fj, c ^/T) /<9T| sat equals zero. Real nonpolar 
solutes display solubility minima at temperatures well 
above the density maximum, largely as a result of at- 
tractive interactions between the solute and water. These 
interactions, not included in the present simulations, can 
be included approximately by assuming they are propor- 
tional to the density of liquid water, as in the van der 
Waals equation-of-state. The resulting chemical p oten- 
tial is — [maIhs — °awPw <|Garde et al\ Jl999). The 
effect of including solute-water interactions on the solu- 
bility is shown in Fig.0 Increasing these interactions sys- 
tematically shifts the maximum in /i^/fcT, out to greater 
temperatures, in agreement with the experimental obser- 
vation of solubility minima at higher temperatures. 

For the hard-sphere solutes the simulation results have 
slightly more curvature at temperatures near the solubil- 
ity minimum than predicted by Eq. (|16[) [Fig. [7] . While 
the fit is accurate, the enhanced curvature suggests the 
heat capacity is not constant as assumed above, but 
is slightly larger at low tempe ratures. Indeed this has 
been observed experimentally l)Gill et all Il985|) . and is 
borne out b y theoretical models of hydrophobic hydra- 
tion as well jAshbaugh et al. L l2002HHummer et ad f2000: 
ISilverstein et all l200l|) . Nevertheless, the temperature 
dependence of the heat capacity is minor, and includ- 
ing it is a complication of secondary importance to the 
interpretations here. We therefore neglect it in our ther- 
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modynamic analysis. 

Fig.|Hlshows how the chemical potential of hard spheres 
in water as a function of temperature with increasing so- 
lute size. In the range of sizes shown, the maximum 
in the chemical potential shifts from temperatures above 
T = 470 K, above the window of temperatures simulated, 
for the 2 A radius solute to lower temperatures with in- 
creasing solute sizes. For cavities not much larger than 
12 A (not shown in figure), the maximum falls below 
T = 260 K, below the simulation window and the nor- 
mal freezing point of water. Thus, for molecular ly sized 
cavities, hydrophobic hydration is opposed by a negative 
entropy, at temperatures below the chemical potential 
maximum, over most of the range of temperatures simu- 
lated. For meso- and macroscopic cavities, however, this 
trend is reversed and hydration is favored by a positive 
dissolution entropy but is insoluble as a result of a dom- 
inating positive enthalpy (demonstrated in the following 
section), dictated by the temperature dependence of the 
surface tension of water. 



radius R — AR . In this case the surface tension is 



D. Surface area dependence of hydration thermodynamic 
properties: curvature corrections and a Tolman length 

In an effort to compare and correlate hydration free 
energies of a variety species, it is common to calculate 
the free energy cost per unit area for hydrating nonpo- 
lar solut e surfaces, also referred to as a m olecular surface 
tension ijHermannL Il977t iTanfordl Il979j) . This molecu- 
lar surface tension, however, is generally not equal to the 
free energy of creating a macroscopic flat interface, in 
part due to curvature and structural differences between 
water at molecular and macroscopic interfaces. Neverthe- 
less, SPT systematically interpolates the surface tension 
between these two length-scal e extremes, and provides 
insight into their relationship l|A shbaugh and P aulaitisl 
l200li iHuang and Chandlel l2000at IHuang et a/J . l200l[) . 

Under the assumption that the pressure contribu- 
tion to the hydration free energy is negligible, an ex- 
cellent assumption for liquid water, a surface tension 
for hydration of a hard-sphere solute is obtained from 
the surface area derivat i ve of the chemical potential 
l|Ashbaugh and Paulaitisl l200l|) . This derivative de- 
pends, however, on the definition of the surface area. A 
natural choice for the solute area is defined by R, and is 
referred to as the solvent accessible surface (SAS) area, 
^4-sas = 47ri? 2 . Differentiating the chemical potential 
with respect to this surface area yields 



7sas(#) 



ox 



d\nR? 



-kT Pw G(R)R 



(17) 



~f(R;AR) 



d 



4tt (R - AR) 2 
1sas(R) 



kT Pw G(R)R 2 
2{R- AR) 



1 - AR/R 



(18) 



More generally, the solute surface can be defined by a 



The SAS tension as a function of R at T = 
300 K is shown in Fig. The surface enthalpy, 
dh%(T )/dA SAS , entropy, ds% (T ) /dA SAS , and heat 
capacity, dc c A * (T ) /dAs AS , are included in this figure. 
For small cavities, all the surface thermodynamic prop- 
erties go to zero as R — > 0. With increasing size, the 
surface tension increases monotonically and approaches 
its asymptotic limit for a flat interface of 7oo = 0.432 
kJ/(mol A 2 ) = 71.7dyne/cm. The other surface prop- 
erties, most notably the heat capacity, approach their 
asymptotic plateaus more slowly with increasing R. Like 
the surface tension, the surface enthalpy monotonically 
increases with increasing solute size. The surface en- 
tropy and heat capacity, on the other hand, vary in dis- 
tinctly different ways for molecular and macroscopic sur- 
faces, indicating changes in the mechanism of hydration 
(ISouthall and Dill [2000) . In particular, the surface en- 
tropy is initially negative beginning from R = 0, consis- 
tent with the experimental thermodynamics of hydropho- 
bic hydration for molecular solutes, reaches a minimum 
at R w 3.3 A and then increases, eventually becoming 
positive as expected from the temperature dependence 
of the liquid-vapor interface. It is curious that the size 
at this minimum is close to the size of the maximum 
of G(R), that is for the most hydrophobic hard sphere 
solute. While the surface heat capacity is positive over 
the entire size range, it reaches a maximum at solute 
radii comparable to the position of the minimum in the 
entropy, suggesting the two are related. Moreover, we 
may infer that the maxima in the surface entropy and 
heat capacity are linked to the breakdown of the aque- 
ous network in the vicinity of the hard-sphere solute as 
observed in simulation studies linking solute- water corre- 
lations to the thermodynamics of hydrophobic hydration 
ijLazaridis and Paulaitisl. Il994l). 

Sharp and coworkers ((Sharp et all Il99lj) suggested 
that rather than relying solely on the SAS to deter- 
mine the molecular surface tension, this tension needs 
to be corrected for the curvature of the molecular inter- 
face to reconcile the difference between molecular and 
macroscopic surface tension. Based on geometric argu- 
ments, they proposed that the radius of a water molecule 
is the length-scale over which this correction must be 
applied. In effect, their work indicates that the van 
der Waals surface, i.e., AR = erww/2 = 1.4 A, pro- 
vides asjop_OTO^_de^mp_tion ofjnole£ukx_TOhite hydra - 
tion l)Jackson and Sternbenl 11994 ISharp et all Il99l|) . 
A schematic illustration of the van der Waals, sol- 
vent accessible, and curvature corrected radii of a hard 
sphere solute in water is given in Fig. 1101 For a 
methane sized solute 7 (R = 3.3A; AR = 0.0A) = 0.300 
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kJ/(mol A 2 ), 30% lower than the macroscopic value, 
while 7(ii = 3.3A;Ai?=1.4A) = 0.521 kJ/(mol A 2 ), 
20% greater than the macroscopic value. While neither 
of these two surfaces gives the macroscopic result, they 
do bracket 7^ suggesting there an optimal intermediate 
value of AR for which the surface tension is size inde- 
pendent. 

Fig. II II shows how the surface tension varies with AR. 
The divergence in 7 (i?; AR) results from the divergence 
in Eq. Q18[) as R — > AR. For solutes larger than 2.5 A, us- 
ing AR — 1 A yields a surface tension that is only weakly 
size dependent. Indeed, 7 (ii = 3.3A; AR = I.OA) = 
0.430 kJ/(mol A 2 ) is in excellent agreement with the 
macroscopic value, suggesting that Honig and coworker's 
geometric estimate of the curvature correction length- 
scale is correct. This argument falls apart, however, when 
we consider the temperature dependence of the Tolman 
length. 

The Tolman length 6, described in Sec. III. Al can be 
accessed by the revised SPT through Eas.l9landlT7l 



7sas (R) ~ 7o 



Substituting this expression into Eq. 
l|Ashbaugh and PauTait 



7(i?;Ai?) ~ 7o 



25/R \ 



1 - AR/R J 



(19) 



flU yields 



(20) 



for the surface tension referenced to a surface displaced 
by AR. Thus, for large R the optimal surface for ob- 
taining a size independent surface free energy is AR = 
28. This Tolman lengt h can be calcul ated from clas- 
sic SPT, see Eq. 10 in (jStillingerL fl973ft . and yields an 
nearly temperature independent 8 w 0.5 A [Fig. I12j. in 
good agreement with the empirical AR at T = 300 K 
obtained above. The revised SPT, however, finds that 
8 is strongly temperature dependent. While the clas- 
sic SPT correctly predicts the magnitude of 8 at low 
temperature, 8 decreases with temp erature and changes 
sign near T = 350 K. Furthermore, l)Moodv and Attardl 
2001) suggest that the Tolman length might change sign 
also for a Lennard- Jones solvent. Thus assuming the 
AR is simply dictated by the size of a water molecule 
leads to fundamentally flawed interpretations of the rela- 
tionsh ip between molecular and macroscopic surface ten - 
sions ijJackson and Sternbergl IT994HSharp et alJ . ll99l[h 
In retrospect, the temperature dependence of the cur- 
vature correction might have been anticipated by the 
entropic differences between hydrating a molecular and 
mesoscopic interface and the significantly different tem- 
perature dependencies of the associated surface thermo- 
dynamic properties [Fig. [5] . The curves observed in this 
figure are simply too rich to be described by a tempera- 
ture independent length-scale. 



E. Entropy convergence and solute size 

When the hydrophobic component of the hydra- 
tion entropies of these molecules is extrapolated to 
high temperatures, the entropies converge to one an- 
other wi t hin a narrow rang e near T = 3 85 K 
llBaldwinl Il98fit iMurnhv et all Il99(t iPrivalovi Il979t 
iPrivalov and Gill Il988|) . This phenomenon of entropy 
convergence is feature of hydrophobic hydration be- 
lieved to be shared between small molecule hydra- 
tion and protein unfolding thermodynamics. Bald- 
win and Privalov noted that this can result from pro- 
portionality of the entropy and heat capacity of hy- 
dration to one another. The most successful ex- 
planations for the convergence temperature for small 
molecules have related the convergenc e temperature to 
the e quation-of-state of pu re wat er llAshbaugh et all 
[ 20021: iGarde arid AshbaughL l200ll iGarde et all Il996t 
iHummer et q/J . ll998j) . IjHuang and Chand ler. 2000_bJ) ar- 
gued that for species larger than R w 10 A entropy con- 
vergence does not occur, and therefore proteins do not 
exhibit this phenomenon . 

In Fig. E| we have plotted the hydration entropies for 
hard-sphere solutes as a function of temperature for so- 
lutes in the size range 2 A < R < 10 A. Indeed, there 
is a broad range of convergence temperatures observed 
at each point where an entropy curve, for a given size, 
crosses that for another solute. For example, the 2 A 
solute entropy intersects the 3 A solute entropy at T = 
410 K, while the 2 A solute entropy intersects the 10 A 
solute curve at T = 300 K, indicating that there is no 
unique convergence temperature. Experimental iden- 
tifications of that entropy convergence, however, have 
largely concerned themselves with solutes similar in size. 
We therefore consider how the convergence temperature 
changes with differential perturbations in the solute size. 
In this case, convergence occurs at the temperature for 
which As A x (R — > R + 8R) = 0. Assuming the hydration 
heat capacity is independent of temperature, T c is deter- 
mined by the relationship 



T c = T exp 



ds%(T ) dcf(T ) 
dAsAs dAsAs 



(21) 



with a size dependence dictated by the relationship be- 
tween ds c A x (To) /dA SA s and dd^ (T ) /dA SA s on the so- 
lute radius in Fig. [5] The differential entropy conver- 
gence temperature curve determined by these expressions 
is shown in Fig. Two points of interest are immedi- 
ately apparent: first the convergence entropy is negative 
and becomes more negative with increasing solute size, 
second the convergence entropy curve more or less forms 
a lower bound on the hydration entropy as a function of 
temperature, although this is approximate. 

One of the implications of Eq. (|21|) is that if the 
entropy is a linear function of the hea t capacity, i.e., 
s ex (Tp) = roc e * (T n ) + 6, as suggested by l|Baldwinl.fl98a 
iMurphv et a71ll99fft . then the convergence temperature 
would be independent of solute size. In Fig. ^] we see 
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that T c has a significant solute size dependence, indi- 
cating that this assumption has limited validity. For 
solutes approaching zero radius, T c plateaus at a max- 
imum. With increasing solute size T c decreases so that 
above R w 8 A it is less than the normal freezing point 
of water. At the intermediate methane radius of 3.3 A, 
however, the convergence temperature is 382K in excel- 
lent agreement with the experimental convergence tem- 
perature of T — 385 K for simple nonpolar gases and 
linear alkanes. 

We note that the convergence temperature as R — > 
plateaus at T = 655 K, above the critical temperature 
of water at 647 K. This unphysical result is due, in part, 
to our extrapolation of the entropies beyond the range, 
260 K to 470 K, to which we fitted Eq. JTSJ. For cav- 
ities small enough that only one water molecule can fit 
inside, the free energy is given by Eq. 10, and entropy 
convergence occurs for the temperature at which 

47T q 

a sa tT c = 1 - — Pw R 3 , (22) 

where a sa t = — (d In pw / ^T) sat is the thermal expan- 
sion coefficient of liquid water along the saturation curve. 
When this criterion is applied, it displays a plateau at 
a physically more realistic temperature of T c « 525 K 
[Fig. El- With increasing solute size, Eq. indicates 
a sudden decrease in T c above R rs 1 A, comparable to 
that obtained on the assumption that the heat capac- 
ity is temperature independent. Above radii of 1.25 A, 
Eq. (|22|l breaks down as multi-particle correlations began 
to play a role in the solute entropy. At this radius, how- 
ever, the convergence temperatures are now within the 
range of temperatures simulated and the application of 
Eq. (|21|l becomes more accurate. It is reasonable then to 
interpolate between the convergence temperatures deter- 
mined by Eq. J2U and Eq. , as indicated in Fig. 

We note that the infor mation theory ijGarde et all 
ll996tlHummer et all\l99^ . with natural simplifying as- 
sumptions, indicates entropy convergence when T c w 
(2a sat ) _1 = 420 K. This corresponds on Fig. [H] to a 
solute radius of R « 2.1 A placing this estimate among 
small-solute theories. Relaxing the assumptions used to 
arrive at this IT criterion lowers the convergence temper- 
ature prediction for methane sized solutes to T 390 K, 
improving agreement with the present result of 382 K. 

The apparent convergence of the entropy change at 
T = 385 K to a value close to zero for a range of hy- 
drophobic solutes and proteins has suggested that the 
hydration of these chemically d i stinct solutes are related 
l|Baldwinl Il986t iMurphv etaR Il990j) . Assuming that 
the temperature dependence for protein unfolding arises 
solely from the exposure of hydrophobic side chains to 
water, phenomenological models have been developed 
which separate out residual temperature independent 
components of the entropy, from contributions such as 
changes in the chain conformation, by extrapolating to 
the convergence temperature. This procedure relies, in 
part, on the assumption that the hydration of surface 



hydrophobic groups is the same in the native and dena- 
tured conformations, and therefore cancels in the entropy 
difference. 

The surprising explanation of the entropy con- 
vergence phenomena that tied it t o the particular 
equat i on of state of liqui d water l| A shbau gh et all 
20021: iGarde and AshbaugH l200lt iGarde et oil Il996t 
Hummer et q?lll998|) resolved an important conundrum 
for our molecular understanding of hydrophobic effects. 
Whether and how this entropy convergence phenomena is 
involved with protein folding is yet an outstanding ques- 
tion. Proteins are complicated molecules participating 
in both hydrophobic and hydrophilic interactions with 
the solution. The widely appreciated point that pro- 
tein folding thermodynamics may be primarily sensitive 
to hydr ation of unfolded configurations is just as im- 
porta nt ijPaulaitis and Pratfl l2002t iPratt and Pohorillel 
2002). Considering unfolded possibilities, the sizes of 
the obvious hydrophobic units are in the range of small 
molecule hydrocarbon solutes. The largest hydro phobic 
side chain — phenyl alanine — is an example. ijPrattt 
2002) emphasized that solution thermodynamic data are 
available for hydrophobic solutes of ju st this size, e.g. for 
benze ne, toluene, and ethyl benzene ijPrivalov and Gill 
119891) . and those data suggest that these solutes exhibit 
conventional entropy convergence behavior. Thus, it is 
a plausible hypothesis that entropy convergence will be 
expressed in protein folding thermodynamics primarily 
through contributions associated with the unfolded con- 
figurations. 

ijHuang and Chandlerl l2000b|) suggested that the hy- 
dration of surface nonpolar groups is better described 
by the hydration entropy of a solute surface compara- 
ble in size to the protein radius on the order of tens 
of angstroms, rather than treating the surface groups 
individually as having sizes comparable to simple hy- 
drophobic units. In this hypothesis it is presumed that 
entropic contributions for hydrating large hydrophobic 
surfaces with attractive dispersion interactions and vic- 
inal polar/charged groups is the same as that for a 
hard repulsive surface. Recent simulations of convex 
methane clusters have found that when realistic attrac- 
tive interactions between water and methane are in- 
cluded, water packs around the cluster methane sites 
j ust as it does around a solitary methane in solution 
ijAshbaueh and Paulaitisl 1200 it IChauL 120031) Moreover, 
Cheng and Rossky ijCheng and Rosskvul998|) found that 
the orientational correlations between water and proxi- 
mal hydrophobic residues on the convex surfaces of the 
bee venom protein, meletin, are comparable to those 
near individual solitary hydrophobic groups in solution. 
These observations suggest that the available configu- 
rational space, and by extension the entropy, for water 
molecules near realistic surface hydrophobic units is the 
same in the folded and unfolded states, supporting the 
assumptions of the phenomenological folding models. We 
note, however, that in the same study Cheng and Rossky 
found that water molecules proximal to hydrophobic 
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residues in flat portions of meletin were more orienta- 
tionally disordered as a result of the difficulties associated 
with maintaining the aqueous hydrog en-bonding; network 
near restrictive solute topologies ijCheng and Rosskvt 
1998). Thus, the applicability of the phenomenologi- 
cal unfolding model may be complicated by the pro- 
tein surface topography and the impact of hydrophobic 
pockets on the overall unfolding entropy. This can in- 
troduce further scatter into m easured folding entropies 
l)Robertson and Murphvl Il997ft . 



IV. SUMMARY AND CONCLUSIONS 

The revised scaled-particle theory bridges the known 
molecular and macroscopic limits by utilizing simula- 
tion information on multi-body correlations in liquid wa- 
ter together with experimental thermodynamic proper- 
ties of pure water to construct a functional form for 
the hydration free energy of hydrophobic hard-sphere 
solutes in water. The classic scaled-particle theory 
llBen-Naim and FriedmMLll967tlMaverlll963tfStillingeil 
11973ft incorrectly predicts an increase in the surface ten- 
sion of water with increasing temperature, and a suspi- 
ciously temperature independent Tolman length which 
does not agree with revised scaled-particle theory obser- 
vations. As a result, application of classic scaled-particle 
theory to hydration free energies is largely a fitting exer- 
cise to obtain the effective van der Waals diameter, ctww, 
of water. Conclusions drawn from this tack have weak 
significance regarding the molecular origins of the hy- 
drophobic effect, and are limited to comparisons of the 
size parameter for water relative to other solvents, ne- 
glecting further molecular detail or specific temperature 
signatures of hydrophobic hydration. 

The revised scaled-particle theory is more success- 
ful, but the success of the scaled particle approach de- 
rives generally from the remarkable fact that the re- 
sults identify a molecular length, near 3.0 A, that pro- 
vides a good joining point for microscopic and macro- 
scopic descriptions. The corresponding results for com- 
parative organic solvents a re less simple (jGrazianol l2003l 
iPratt and PohorilleL ir9921. That micro-macro joining ra- 
dius exhibits interesting temperature variation; an accu- 
rate description of those temperature variations is an im- 
portant part of the higher fidelity of the revised scaled- 
particle results. The revised scaled-particle theory re- 
produces the solubility minimum behavior for small hy- 
drophobic solutes, and demonstrates significant changes 
in the hydration mechanism of hard-sphere solutes with 
increasing solute size. Specifically, the hydration thermo- 
dynamics of small solutes is predominantly entropic at 
room temperature, but the hydration of mesoscopic cav- 
ities is entropically favorable and opposed by a dominat- 
ing hydration enthalpy. While it is tempting to describe 
these changes in hydration thermodynamics in terms of 
aqueous hydrogen-bonding near the hydrophobic entity 
— and that can be plausible in the appropriate theoret- 



ical setting — the scaled-particle theory provides little 
in the way of information on the integrity of hydrogen 
bonded networks. 

Nevertheless, the revised scaled-particle theory does 
provide thermodynamic information that challenges phc- 
nomenological views of hydrophobic effects, particularly 
the cherished iceberg hypothesis. Whereas the iceberg 
hypothesis suggests that local freezing of water molecules 
in the vicinity of hydrophobic solutes is a source for the 
negative hydration enthalpies, we find that at room tem- 
perature the hydration of solutes comparable in size to 
simple nonpolar gases is actually unfavorable from an 
cnthalpic as well as an entropic standpoint. The exper- 
imentally determined favorable enthalpies of solution of 
hydrophobic species then are a consequence of attractive 
solute-water interactions and not enhanced water-water 
structuring. 

On a molecular level there is a surface that maps 
macroscopic surface tensions to molecular values. This 
reduces the reconciliation of molecular and macroscopic 
values of the surface tension to a program of finding the 
appropriate dividing surface. The utility of that pro- 
gram rests on the optimistic expectation that the Tolman 
length is largely temperature independent. But that the 
Tolman length was found to have a significant tempera- 
ture dependence in water, changing from positive to neg- 
ative at T ss 350 K, a possibility anticipated by Stillinger 
llStillingeilll97l . As a result then, though the optimal 
surface for the description of hydration may be approx- 
imated by the solute molecular surface at low temper- 
ature s ijAshbaugh et all Il999t lAshbaugh and Paulaitisl 
l200lj) . with increasing temperature this optimized sur- 
face moves out to the solvent accessible surface at T w 
350 K, and ultimately extends beyond this surface at even 
higher temperatures as a result of the nontrivial temper- 
ature dependence of the hydration thermodynamics for 
molecular-sized solutes. 

Finally, the revised scaled-particle theory provides de- 
tailed information on the curious entropy convergence 
behavior observed for small molecule solutes and the size 
dependence of the convergence temperature. A suitably 
defined entropy convergence temperature retreats below 
the freezing temperature of water for hard spheres the 
size of globular soluble proteins. But heterogeneity of 
protein-water interactions and of sizes of hydrophobic 
units also contribute importantly to experimental blur- 
ring of entropy convergence behavior in protein unfold- 
ing thermodynamic data. Equally important, entropy 
convergence behavior for protein folding thermodynam- 
ics may be primarily expressed through contributions as- 
sociated with the unfolded configurations, and due to hy- 
dration of hydrophobic side chains of size corresponding 
to studied small molecule solutes. 
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FIG. 1 Product of the liquid-vapor surface tension, 7, and 
the the isothermal compressibility, k,t = —y (lip") > f° r sev " 
eral liquid solvents at low temperatures so that the density 
of the coexisting vapor is low. Even though the individual 
factors differ substantially in magnitude, this product is a 
length characteristic of the liquid, accessible on the basis of 
macroscopic measurements, and can be taken as proportional 
to a molecular correlation length. The interesting observa- 
tion here is that temperature dependence of this correlation 
length is qualitatively different for liquid water than for these 
organic solvents. 



FIG. 2 The box on the left shows a configuration of water 
molecules taken from a simulation of liquid water at the den- 
sity of the liquid in coexistence with vapor at 300 K. The 
box on the right shows hard spheres of diameter 2.8 A that 
can be successfully placed into the configuration on the right 
without overlap of the van der Waals volume of the water 
molecules. The insertion probability po is determined as the 
volume accessible to centers of the purple spheres divided by 
the geometric volume of the box. 
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FIG. 3 Cavity contact function for water at T = 300 K at 
the liquid saturation conditions. The points are obtained by 
differentiation of the simulation cavity insertion probabilities. 
The dashed lines are obtained from Reiss's original SPT the- 
ory predictions for a hard sphere solvent, Eq. 110L using ef- 
fective water hard sphere diameters between oww = 2.6 A to 
3.0 A in 0.1 A increments. The solid line is obtained by differ- 
entiation of the revised SPT, Eq. I15L fitted to the simulation 
results between i? s im = 2.5 A and i? m acro = 3.5 A. 
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FIG. 4 Cavity-water oxygen radial distribution function for a 
3 A cavity at T = 300 K. The thin solid line indicates the ra- 
dial distribution function, while the dashed line indicates the 
radial integral N(r) = JT pwg(A)47rA 2 dA. The occupation 
of the first hydration shell, corresponding to the first mini- 
mum in g(r) at 5.1 A, is 17.8 water molecules as indicated 
by the thick horizontal line. Note that the first minimum, 
which physically discriminates between first and succeeding 
hydration shel ls, is mild, and structurin g of outer hydration 
shells is weak llPratt and Pohorillel . ll993Tl . These features are 
in good qualitative agr eement with the predictio ns of the 
Pratt-Chandler theory iPratt and ChandleU Il977f) . though 
that specific theory has been substantially amended iPrattl 
2002). 
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FIG. 5 Cavity contact function for water as a function of tem- 
perature along the liquid saturation curve determined by the 
revised SPT, Eq. JUJ, with _R aim = 2.5 A and R macfo = 3.5 A. 
Results are shown between T = 260 K to 460 K in 20 K incre- 
ments. Notice that the length defined by the maximum of this 
curve decreases with decreasing density following increasing 
temperature along the saturation curve. This is consistent 
with the suggestion of Fig. 
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FIG. 6 Excess chemical potential, enthalpy, and entropy of 
a methane sized hard sphere solute (R = 3.3 A) in water as 
a function of temperature along the saturation curve. The 
points are the explicit simulation results for the chemical po- 
tential. Error bars are comparable in size to the points. The 
curves for the excess chemical potential, enthalpy, and en- 
tropy are labeled in the figure. The curves were determined 
under the assumption that the heat capacity is independent 
of temperature [Eqs. 1161 1. 
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FIG. 7 Solubility of a 3.3 A hard-core + weak attractions 
solute as a function of temperature with increasing strength 
of attractive interactions modeled as in the van der Waals 
equation-of-state, pi™ = [Ma'Ihs — a swpw- The points are 
simulation results from particle insertion probabilities. The 
solid curves are the solubilities with lower curves indicating 
increasing attractive interactions, asw > 0. The dashed curve 
indicates the maxima in n^/kT with increasing interactions, 
corresponding to minima in the Ostwald solubility, defined as 
K eq = exp [-[if (aq) /kT] [Eq. @]. 
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FIG. 8 Excess solute chemical potential as a function of tem- 
perature for solutes of varying size. The solid lines corre- 
spond to the revised SPT model. The solid circles, open cir- 
cles, and crosses correspond to explicit molecular simulation 
results simulation results for the 2 A, 3 A, and 4 A radius 
solutes, respectively. Estimated statistical errors are smaller 
than the plotting symbols. The dashed line indicates the locus 
of chemical potential maxima, where = 0, with changing 
cavity size. 
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FIG. 9 SAS area derivatives of the hard-sphere solute hy- 
dration thermodynamics as a function of solute radius at 
T — 300 K. The thin solid, long-short dashed, and short 
dashed, and long dashed lines correspond to dfi c \ /OAs as, 
dhT/OAsAS, TdsT/dA S As, TdcT/dAsAS, respectively. The 
thick horizontal line indicates the macroscopic surface tension 
for a flat surface. 



FIG. 11 Surface tension for hydration of hard-sphere solutes 
as a function of solute radius at T = 300 K employing dif- 
ferent definitions of the solute surface. The thick solid line 
corresponds to the surface tension determined by the deriva- 
tive with respect to the solvent accessible surface area defined 
by the radius R [Eq. UlTI l . The lines above the baseline SAS 
surface tension indicate the effect of increasing AR in 0.25 A 
increments from 0.25 A to 1.25 A [Eq. JUJ]. The thick dashed 
line corresponds to AR = 1 A. 
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FIG. 10 Alternative definitions of the radius of a hard sphere 
cavity. The solvent accessible radius, R, is given by the dis- 
tance of closest approach between the center of the cavity and 
the water oxygen center. The radius R — oww/2 demarks the 
van der Waals boundary of the cavity at ctaa/2. R — AR. 
locates a neighboring surface that might provide a curvature 
corrected surface tension. 
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FIG. 12 The curvature correction, S, as a function of temper- 
ature along the saturation curve of water. The points corre- 
spond to the values determined by the fit of Eq. (11411 to the 
simulation free energies. The solid line is a guide to the eye for 
the fitted results. The dashed line corresponds to the classic 
SPT prediction for the Tolman length, Eq. 10 in (IStillingeit 

inn 



R=10A 9A 8A 7A 





0- 


2A 
3A 


i 


-100 - 


4A 


"o 






E 


-200 


5A 








s < 

C-3 


-300 


6A 




FIG. 13 Entropy of hydrophobic hydration as a function of 
temperature for solutes in the size range 2 A< R <10 A m 
1 A increments. The dashed lines are the excess entropies of 
hydration. The open circles are the convergence temperatures 
for consecutive solutes, i.e., s A x (R) = s\(R + lA). The thick 
solid line indicates the entropy convergence temperature in 
the limit of infinitesimal perturbations in R. 
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FIG. 14 Variation of the entropy convergence temperature 
with increasing hard-sphere radius. The thin solid line is the 
convergence temperature determined under the assumption 
the heat capacity is independent of temperature. The thick 
solid line is the exact entropy convergence temperature for 
R < aww/2 from Eq. I|22|l , The dashed line smoothly inter- 
polates between the exact and constant heat capacity curves 
at 1.25 A and 3.3 A, respectively. The filled circle indicates 
the entropy convergence temperature of a methane-sized so- 
lute (T c = 382 K) . The open circle indicates the entropy con- 
vergence temperature based on the IT criterion (T c = 420 K) . 



